Transient growth of Kelvin waves on quantized 

cn 1 vortices 

O ' 
O 

CN ■ C. F. Barenghi 1 , M. Tsubota 2 , A. Mitani 2 and T. Araki 2 

C3 ' 
33 



(N 



1 School of Mathematics and Statistics, University of Newcastle, Newcastle, NE1 7RU, UK 
2 Physics Department, Faculty of Science, Osaka City University, Japan 



Introduction 



[ Quantized vortex lines in helium II can be destabilized by a sufficiently large 

■ normal fluid velocity which is parallel to the vortex lines (Donnelly - Glaber- 

l/"") . son instability). We study what happens if the driving normal fluid is not 

steady but oscillates periodically with time. We show that in certain situa- 
^ ■ tions, although individual perturbations decay for t — » oo ; a state of transient 

growth of Kelvin waves is possible, which may explain observations of a small 
+^ \ and not reproducible longitudinal component of mutual friction. 

a 

C3 
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q , Recent work on the hydrodynamics of helium II has highlighted the 

importance of the Donnelly - Glaberson instability. The instability takes 
the form of Kelvin waves (rotating helical displacements of the vortex core) 
which grow exponentially along quantized vortex lines. If this happens, 
the length of the vortex lines increases and energy is transferred from the 
normal fluid to the superfluicP. The Donnelly - Glaberson instability was 
discovered^ when a small heat current was applied in the direction parallel 
to quantized vortices in rotating helium II. Figure 1 shows this rotating 
counterflow configuration. A channel is closed at one end and open to the 
helium bath at the opposite end. At the closed end a resistor R dissipates a 
known heat flux Q. The heat is carried away from the resistor by the normal 
fluid whilst the superfluid flows in the opposite direction to assure that the 
net mass flux is zero. In this way a counterflow velocity V is set up which 
is proportional to the applied heat flux. If the channel is rotated about 
its axis at angular velocity O, in the absence of heat flux there would be a 
uniform array of superfluid vortices of areal density L = 2Q/r aligned along 
the rotation axis where T is the quantum of circulation. In the presence 
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of a heat flux, the normal fluid velocity v n is parallel to the vortices and, 
if Q is large enough, the Donnelly - Glaberson instability takes place and 
Kelvin waves grow exponentially with time. When the amplitude of the 
Kelvin waves becomes of the order of the intervortex spacing b ~ L -1 / 2 , 
vortex reconnections take place and, according to a recent calculation^, the 
uniform array of vortices may turn into a state of polarized turbulence. 

It must be stressed that, although the rotating heat transfer configu- 
ration of Figure 1 is the best to study the Donnelly - Glaberson instability 
in a controlled way, there are other experiments in which the Donnelly - 
Glaberson instability plays a key role. For example, if helium II is made 
turbulent using a towed gricP or rotating propellers®, both normal fluid and 
superfluid are turbulent (the former consisting of eddies of various sizes and 
strengths, the latter consisting of a disordered tangle of quantized vortices). 
In this case the instability may take place throughout the turbulent flow, 
provided that the component of v n locally parallel to a superfluid vortex 
filament exceeds a critical value 

The aim of this paper is to remark that a similar Donnelly - Glaberson 
instability can be excited if v n is not steady but varies harmonically with 
time. The natural way to do this is to modify the configuration of Figure 1 
and connect the resistor R to an AC (rather than a DC) power supply. For 
the sake of brevity we shall refer to the old and new configurations as the 
"DC instability" and the "AC instability" respectively. 

The plan of the paper is the following: after deriving the DC instability 
in Section 2, we describe two versions of the new instability in Sections 3 and 
4. The results are discussed in Section 5 in view of previous experiments. 



2. The DC instability 

Let us consider a quantized vortex which points in the z direction and 
assume that its shape is perturbed in the form of a small amplitude Kelvin 
wave of wavelength A. A point on the vortex has Cartesian components 

s(£) = [A cos (0); A sin (<£);£] (1) 

where 4> = k£ — cut, t is time, and u, A and k = 2tt/X are respectively the 
angular frequency, the amplitude and the wavenumber of the wave. Assum- 
ing Ak « 1, it is easy to verify that the parameter £ is the arc length. In 
fact s' = ds/dt; = [—kAsm(<f));kAcos((ft);l] is the tangent to the curve at 
the point s and the length A measured along the curve is 

A= /^Wls'te^-s'te,)! = / cWfc^ + lwf. (2) 
Jo Jo 
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In the absence of friction and using the local induction approximation^, 
the self-induced velocity of the vortex at the point s is 

— = v se// = i/ s s' x s", (3) 

where u s = (T/(4n)) log (b/a) and a is the vortex core radius. From Q and 
(|3"|) we obtain the dispersion relation of Kelvin waves 

u) w v s k 2 . (4) 

In the presence of friction the equation of motion is ^ 
ds 

-jj; = Vself +« s ' x ( v n - V s - V 8e if) (5) 

where a is a known mutual friction coefficient and v n and v s are imposed 
normal fluid and superfluid velocity fields. Choosing the z axis along the 
counterflow channel as in Figure 1, we have v n — v s = [0; 0; V] where V is the 
applied counterflow velocity. Allowing the amplitude of the wave to change 
with time and using (jlj) we obtain 

A A 

— = a(kV - u s k 2 )A, (6) 

If V is constant the solution of equation © is 

A(t) = A(0)exp(pt), (7) 

where A(0) is the initial infinitesimal amplitude of the wave due to thermal 
or mechanical vibrations and a = a(kV — v s k 2 ) is the growth rate. The 
instability occurs if a > 0. Therefore a Kelvin wave of wavenumber k is 
unstable in the region < k < ko where ko = V/v s and the maximum 
growth rate a c = aV 2 /{Av s ) occurs at wavenumber k c = V/{2v s ) = ko/2. 
The temperature of the helium bath determines a and the applied heat 
current determines V, so the Kelvin wave which grows exponentially faster 
and dominates all the other waves has wavelength A c = 47r^ s /y. 



3. The AC instability: second sound 

In second sound the quantity V reverses sign during a cycle. This can be 
achieved by using Nuclepore transducers, which are oscillating membranes 
covered by holes which are so small (sub-micron size) that the viscous normal 
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fluid cannot move through them. In this way an oscillating counterflow 
velocity is set up and Eq. © becomes 

dA 

— = a (kV cos{u t) - v s k 2 )A, (8) 

Letting t' = u^t and \x = u s k/Vo the solution of Eq. (jHJ) is 

A{t') = A(0)efV\ (9) 

where 

/(0 = — [sm(t')-/^], (10) 

The amplitude A of the perturbation grows if / > 0, but, since | sin(t') | < 
1, we have / < for sufficiently long times. We conclude that infinitesimal 
perturbations of any wavenumber k decay exponentially for t — > oo. How- 
ever, although there are no instabilities for t — > oo, for a limited initial time 
helical waves can actually increase their amplitude. This "transient growth" 
takes place in the time interval such that sin(i') > fit' that is to say for 

v s kujQ 

The effect can be significant (it can last for an appreciable time) at low 
second sound frequency (ljq « 1), high second sound amplitude (Vb >> 1) 
and long wavelength (k « 1), because v s w 10 _3 cm 2 /sec << 1 (taking 
ao ~ 10 _8 cm and O w lrad/sec). Note that the longest wavelength depends 
on the size of the channel. The important point to notice is that a great range 
of wavenumbers can be affected by transient growth. Since mechanical and 
thermal perturbations are created all the time, if these perturbations grow 
before decaying to zero at t — > oo, then each vortex line will have a certain 
amount of small amplitude Kelvin waves, depending on the particular set 
up of the experiment (ujo, Vq and the size of the channel) By attenuating 
the second sound which propagates in the direction along the vortices, these 
Kelvin waves may lead to think that there exists a longitudinal coefficient 
of mutual friction. 



4. The AC instability: oscillating counterflow 



Perhaps the simplest way to produce a time dependent heat flux Q is 
to apply an alternating voltage to the resistor R which generates the heat 
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current. Since Q is proportional to the square of the applied voltage, Eq. (JSJ) 
must be replaced by 

dA 

— = a [kV cos 2 (uj t) - u s k 2 ]A, (12) 
Preceding as before we have 

[cos 2 (0 - fi]A, (13) 



dA akVo 2 



and the solution is 
where 



dt' ujq 

A{t') =A(0)eW\ (14) 

/(*') = ^[(l-2M)t' + ^sin(2t')], (15) 

The perturbations grow if fit') > 0. The long term behaviour of f(t') 
depends on the first term in Eq. (|15|) which increases with t', not on the 
second which is always less than unity. If the first term is positive (1— 2/i > 0) 
then A grows for t — > oo; this happens for /i < 1/2, that is for < k < 
which is the same condition as for the DC instability with V$ replaced by 
Vq/2. Transient growth is possible for a limited time even if 1 — 2fj, < 
provided that is for t < l/(2u;o|l — 2//|). 



5. Conclusion 

We have shown that there is an "AC" version of the "DC" Donnelly - 
Glaberson instability, and that can it be realized in a controlled way in the 
laboratory, either using second sound Nuclepore transducers or an oscillating 
counterflow. We have determined the boundaries between the stable and the 
unstable regions and found that in the linearly stable regions, although each 
individual perturbation decays as t — > oo, a state of transient growth is 
possible which depends on the frequency and the amplitude of the second 
sound wave and on the length of the channel. This state of transient growth 
of Kelvin waves would produce attenuation of second sound in the direction 
parallel to the quantised vortices and would be interpreted as evidence for a 
longitudinal component of mutual friction. The possibility of the existence 
of a longitudinal friction coefficient (called B", B z or 7) was considered in 
the earlier helium II literature^. Experiments^ concluded that B", if it 
exists at all, is very small (at the most 500 times smaller than B), is not 
reproducible and may depend on the amplitude of the second sound wave 
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Fig. 1. Rotating counterflow channel configuration. 



used to measure it, which is consistent with our finding. We plan to continue 
this investigation using numerical simulations. 

C.F.B. is grateful to the Royal Society for supporting this study. 



REFERENCES 

1. D.C. Samuels and D. Kivotidcs, Phys. Rev. Lett. 83, 5306 (1999). 

2. D. K. Cheng, M.W. Cromar and R.J. Donnelly, Phys. Rev. Lett. 31, 433 (1973). 

3. R.M. Ostermeier and W.I. Glaberson, J. Low Temp. Phys. 21, 191 (1975); 

4. M. Tsubota, T. Araki and C.F. Barenghi, Phys. Rev. Lett. 90, 205301 (2003). 

5. L. Skrbek, J.J. Niemela and R.J. Donnelly, Phys. Rev. Lett. 85, 2973 (2000). 

6. J. Maurer and P. Tabeling, Europhys. Lett. 43, 29 (1998). 

7. C.F. Barenghi, D.C. Samuels, G.H. Bauer and R.J. Donnelly, Phys. Fluids 9, 
2631 (1997). 

8. C.F. Barenghi, S. Hulton and D.C. Samuels, Phys. Rev. Lett. 89 275301 (2002). 

9. P.G. Saffman Vortex Dynamics, Cambridge University Press, (1992). 

10. K.W. Schwarz, Phys. Rev. B 38, 2398 (1988). 

11. I.M. Khalatnikov, An Introduction To The Theory Of Superfluidity, W.A. Ben- 
jamin Inc. (1965). 

12. H.A. Snyder, Phys. Fluids 6, 755 (1963). 



